Priedel oscillations in disordered quantum wires: Influence of e-e interactions on the 

localization length 
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The Friedel oscillations caused due to an impurity located at one edge of a disordered interacting 
quantum wire are calculated numerically. The electron density in the system's ground state is 
determined using the DMRG method, and the Friedel oscillations data is extracted using the density 
difference between the case in which the wire is coupled to an impurity and the case where the 
impurity is uncoupled. We show that the power law decay of the oscillations occurring for an 
interacting clean ID samples described by Luttinger liquid theory, is multiplied by an exponential 
decay term due to the disorder. Scaling of the average Friedel oscillations by this exponential 
term collapses the disordered samples data on the clean results. We show that the length scale 
governing the exponential decay may be associated with the Anderson localization length and thus 
be used as a convenient way to determine the dependence of the localization length on disorder 
and interactions. The localization length decreases as a function of the interaction strength, in 
accordance with previous predictions. 

PACS numbers: 71.55.Jv, 71.55.-i,73.21.Hb 



I. INTRODUCTION 

The interplay between repulsive interactions and dis- 
order in low dimensional systems, and their influence on 
the conductivity, were the subjects of many studies in 
recent years. Some of this interest was motivated by the 
experimental observations of a crossover in the temper- 
ature dependence of the conductance of low density two 
dimensional electrons from an insulating like dependence 
at low densities to a metallic one at higher densities-. 
Nowadays it is generally accepted that even if such a 2D 
metal insulator transition exists, it must be related to the 
spin degree of freedom^ and therefore absent for spinless 
electrons. 

It seems therefore clear that for spinless one- 
dimensional systems no metal insulator transition is ex- 
pected for repulsive interactions, although for a certain 
range of attractive interactions a delocalized regime was 
found in several studiesS.. Nevertheless, it was shown that 
there might be a certain strong disorder and interaction 
regime, in which the localization length, or other prop- 
erties usually related to it such as the persistent current, 
increase^. A sample dependent increase in the localiza- 
tion length was also reported for weaker values of disor- 
der and repulsive interactions for longer (of order of 100 
sites) wires^. 

On the other hand, several analytic studies^ have con- 
cluded that the localization length decreases monotoni- 
cally with increasing repulsive interaction. Using either 
renormalization group^ or self consistent Hartree Fock^ 
methods it was shown that the localization length, renor- 
malized by the interaction, scales as 



- (Co) 



l/(3-2g) 



(1) 



where is the localization length of the free electron 
system, and g is the TLL (Tomonaga-Luttinger liquid) 
interactions parameter with g = 1 for non-interacting 



electrons. Since for repulsive interactions g decreases as 
a function of the interaction strength, one finds that the 
localization length always decreases as a function of the 
interaction strength. 

One must be careful though to differentiate between 
weak and strong interaction strength. A clean one dimen- 
sional system of spinless fermions on a lattice undergoes 
a metal-insulator phase transition between a TLL and a 
charge density wave (CDW) as a function of the interac- 
tion strength. This transition, caused by umklapp pro- 
cesses, is exhibited for commensurate fillings. Once dis- 
order is turned on, the TLL transport properties change 
drastically. For more than a decade it is well knownS that 
the conductivity of a TLL wire vanishes in the presence 
of impurities, thus a metal- insulator transition as a func- 
tion of interaction strength no longer exists. Yet, there 
is a difference between the two phases, since the TLL 
is replaced by an Anderson insulator, while the CDW 
phase may remain a Mott type insulator, or become an 
Anderson insulator— lii. 

In this paper we investigate numerically the regime of 
the Anderson insulator caused by adding disorder to the 
TLL phase. We study the effect of the interplay of weak 
interactions and disorder on the behavior of the Friedel 
oscillations in a wire due to its coupling to an impurity 
at its edge. Strictly speaking we probe the exponential 
decay of the Friedel oscillations as a function of disor- 
der and interactions, but for weak disorder this decay 
length is equivalent to the localization length. It is im- 
portant to note that the extraction of the localization 
length for interacting systems is plagued with difhcul- 
ties. The straightforward method of measuring the decay 
length of the envelope of the single electron state has no 
direct translation to a many electron state. Nevertheless, 
one would prefer to stick to a ground state property of 
the system, since the calculation of excited state depen- 
dent properties such as the conductance is computation- 
ally taxing. The sensitivity to boundary conditions (i.e.. 
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persistent current) which is the natural candidate for a 
ground state property is problematic since it incorporates 
both interaction corrections to the localization length as 
well as interaction corrections to the inverse compress- 
ibility of the system^^. Thus, the study of the influence 
of interaction on the Friedel oscillation in the Anderson 
phase is not only interesting on its own account, but it 
establishes a new numerical method using a ground state 
property which is convenient for a direct evaluation of 
the localization length for not too strong disorder. Using 
this method we show that the localization length as a 
function of the interaction strength decreases, in corre- 
spondence to Eq. 

In general, the study of a dot (or impurity) coupled 
to a one-dimensional lead, has been shown to shed some 
light over the physics of the lead. Certain thermody- 
namic observables, such as the occupation of the impu- 
rity leveli^ii^ii^ and the corresponding electron density 
changes in the leadi^iii, were recently used to analyze 
different wire properties, such as the strength and form 
of the interactions, and even the wire's phase (e.g., TLL 
vs. CDW). In a similar fashion, we show how the electron 
density of a disordered wire, coupled to an impurity level, 
can be used in order to extract its localization length. 

Once a single-level impurity (dot) is coupled to a clean 
metallic system the density of electrons in its vicinity 
oscillates with a 2kp period, and the envelope of the os- 
cillations decays as a power law of r, the distance from 
the impurityl^. For non-interacting systems the pertur- 
bation of the density in the vicinity of the impurity de- 
pends on the dimensionality, d, of the system, and can 
be expressed as 



Sp{r) = A 



cos{2kpr + rj) 



(2) 



where the coefficient A and the phase shift rj do not de- 
pend on r. These oscillations are the famous Friedel Os- 
cillations (FO), which have been observed experimentally 
during the last decade using various techniques, such as 
scanning tunneling microscopy in low temperatureslS' and 
X-Ray diffraction^-. 

Whereas for higher dimensions {d > 2) Eq. is in 
general true even in the presence of interactions, this is 
not the case for ID systems. For the TLL phase, using 
field theoretical approaches, it was shown™ that the 
dependence is replaced by a different power law, . For 
the non- interacting case g = 1, it leads to the expected 
decay, while for repulsive interactions g < 1 and thus 
a slower decay of the FO envelope is expected. 

From Eq. it is clear that the observation of the den- 
sity fluctuations, either experimentally or numerically, 
is easier at short distances in the vicinity of the impu- 
rity. When disorder is also introduced, this distance be- 
comes even shorter since there are also density fluctua- 
tions caused by the disorder. Yet, in common experimen- 
tal ID situations disorder is usually present. Therefore, 
although the presence of disorder hampers observing the 
FO, it is beneficial to develop a method to tease the FO 



out of the density fluctuations of a disordered system. 

The paper is organized as follows. In the following 
section we present the system's many particle Hamilto- 
nian and the diagonalization method. We also describe 
a simple method used to extract the FO data of disor- 
dered samples. In the current paper we restrict ourselves 
to the weak interactions regime (TLL), and the results 
are presented in section 3. Results for the CDW phase 
(strong interactions), which show quite different physics, 
will be presented elsewhere^. In the last section we dis- 
cuss the results mostly by a qualitative comparison to 
previous predictions, and offer some possible experimen- 
tal realizations. 



II. METHOD 

A. Hamiltonian and diagonalization method 

The system under investigation is composed of a spin- 
less one dimensional electrons coupled to an impurity in 
one end. We model the one dimensional wire by a lat- 
tice of size L with repulsive nearest neighbor (NN) in- 
teractions and with an on-site disorder. The system's 
Hamiltonian is thus given by 

L L-1 
H^^re = ^ ejcjcj - t ^ (c]Cj + i -f- h.C.) (3) 

L-1 



where ej are the random on-site energies, taken from a 
uniform distribution in the range [—14^/2, 14^/2], / is the 
NN interaction strength (/ > 0), and t is the hopping 
matrix element between NN, henceforth taken as unity. 
c| {cj ) is the creation (annihilation) operator of a spinless 
electron at site j in the wire, and a positive background 
is included in the interaction term. 

Without the disorder term, a similar system - in the 
limit L ^ oo and with periodic boundary conditions - 
has a well known exact solution. Depending on the in- 
teraction strength, the wire can be either metallic or in- 
sulating. The metallic phase is described by TLL, oc- 
curring for / < 2t, and the insulating phase, in which 
/ > 2t, is a CDW. Previous studies which used wires of 
the order of a few hundreds sites have shown a similar 
phase diagrami^ii^ even when employing open boundary 
conditions. In order to stay in the TLL (Anderson insu- 
lator) regime for the clean (disordered) case, we restrict 
ourselves to the range < I < 2t. 

Introducing an impurity at one end of the wire results 
in adding the following term to the Hamiltonian: 



imp — ^O^^o'^O 



y(cSci - 

\){c\ii 



h.C.) 



(4) 
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where eg describes the impurity strength, and V is the 
hopping matrix element between the impurity and the 
first lead site. Along this paper we use cq ^ W and 
V = t. 

The resulting Hamiltonian H — Hwire + tiimp de- 
scribes a disordered one dimensional wire of length L 
{1 < j < L), which is coupled to a single level at one of 
its edges (j — 0). Practically the j = site is equivalent 
to any other site, except for having a constant onsite en- 
ergy, whereas the other sites have energies drawn from 
a distribution with a zero average over different realiza- 
tions. 

The Hamiltonian H was diagonalized using a finite-size 
DMRG methodi^i^^, and the occupation of the lattice 
sites were calculated, for different values of eg, W and /. 
The size of the wire was up to i = 500 sites. During 
the renormalization process the number of particles in 
the system is not fixed, so that the results describe the 
experimentally realizable situation of a finite section of 
a ID wire which is coupled to a dot and to an external 
electron reservoir. 



B. Extracting the Friedel oscillations decay 

When no disorder is present {W = 0), Hwire has a 
particle-hole symmetry, and the particle density of the 
wire's ground state is flat, with filling factor n — 1/2. 
In this case 2kF = tt and the oscillating part of Eq. 

wire+imp 



alternates according to (— I)-'- Denoting by n 
f^^mre-j ^]^g electron density at site j of the wire when 
coupled (not coupled) to the dot, one has n™''^ — n — 
1/2 for any j. Clearly this is not the case in the presence 
of the impurity, and the effect of the impurity is measured 
by Nj = — n. A typical result of Nj, showing 

the 2kF oscillations caused by the impurity at j — 0, is 
presented in Fig. [TJ 
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FIG. 1: A typical form of tlie FO (eo = 10 and L = 280) 
without disorder and witliout interactions. The impurity is 
located at j = and the population of the first 20 lead sites 
is shown. 



When W 7^ 0, on the other hand, although the aver- 
age filling factor is still n ~ 1/2, there is no local sym- 
metry between particles and holes, and the disorder ef- 
fects are seen in the fluctuations of the electron density. 
The density oscillations generated by the additional im- 
purity are then difficult to discern, since in a distance of 
a few lattice sites from the impurity the disorder fluctua- 
tions are dominant. A typical result of Nj together with 
Nj = nj^'^'^ — n (the electron density of the disordered 
wire without an impurity), is shown in the upper panel 

of Fig. m 
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FIG. 2: A typical FO for a disordered sample with L = 280, 
VF = 0.1 and eo = 10 (without interactions). The upper 
panel shows Nj (circles) and A'^' (squares), and the lower 
panel presents the difference between them {ANj). The FO 
are observed much better using AA'j instead of Nj. 

However, the influence of the impurity can be observed 
by isolating the density fluctuations created by the disor- 
der. This is achieved by comparing the electron density 
of the two cases shown in the upper panel of Fig. [21 i.e. 
one with the additional impurity and the other without 
it, for every disorder realization. Averaging over realiza- 
tions is thus done for 



AA^,- = N. - N. 







(5) 



instead of just averaging over Nj. The curve of ANj in 
the lower panel of Fig. [2] is for the same realization as in 
the upper panel. It is obvious that the FO which were 
hardly seen for Nj become clear once ANj is considered. 



III. RESULTS 

We begin with the results for a clean sample (i.e., W = 
0), so that the density of the lead (disconnected from the 
impurity) is half everywhere. As was mentioned above, 
for a non-interacting wire (/ = in our model) , the decay 
of the oscillations is inversely proportional to the distance 
from the impurity, and one expects to get, except for very 



short distances from the impurity, ANj 



N, 
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A{—l)^j ^, where the amplitude A does not depend on 
J- 

However, the ampHtude A does depend on eg, as pre- 
sented in Fig. [21 The hmits of eo ^ and sq —^ oo are 
weh understood, because in both of them the impurity 
does not play any role, the lead has a hard wall boundary, 
and the particle-hole symmetry imposes that the FO am- 
plitude goes to zero. For finite values of eo, the behavior 
of the amplitude is found to obey the relation 

for which a complete proof is given in the appendix. The 
correspondence between the numerical results and this 
formula, as shown in the inset of Fig. [31 is excellent. 
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FIG. 3: The FO decay (in log-log scale) for different values 
of eo (shown next to the curves) for clean samples without 
interaction. The slope - representing the decay exponent - is 
constant, and the only effect of eo is a change of the amplitude 
A. The curves for eo = 2.0, 10.0 are drawn with symbols. 
Inset: the dependence of A on eo for V = t together with the 
exact formula Eq. (|6]) which is derived in the appendix. 

We now move to the interacting case. For < I < 2t, 
i.e. when the fermions in the lattice are described by the 
TLL theory, the decay is expected to be proportional^^ 
to j^^ ■ In our model the TLL parameter g is given by^ 



2cos-i[-^/(2t)]' 

For non-interacting particles one gets g = I so that the 
oscillations decay as while in the interacting regime 
a monotonic decrease of g toward the limit g = 1/2 oc- 
curs as a function of interaction strength. Thus, as / 
becomes stronger, g decreases, and a slower decay is pre- 
dicted. This trend is seen in the DMRG results presented 
in Fig. H 

In the inset of Fig. [H the results obtained for g by fit- 
ting the FO decay of a 500 sites wire, to the predicted 
decay of , are presented together with the theory pre- 
diction for g{I) of Eq. ([7]). As can be seen, the results 
are in good accordance with the theory for interaction 
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FIG. 4: The DMRG results for the FO decay in log-log scale 
for I = 0,0.25,0.75 and 1.5 (bottom to top), with eo = 10 
and L = 280 and without disorder. As / increases, the decay 
gets slower. Inset: The interaction parameter g as found by 
fitting the FO decay to (symbols), together with the the- 
ory prediction Eq. (|7]l (line). These results were obtained by 
taking L = 500 and eo = 1. 

strength I/t ^ 1. Similar results were obtained using 
other implementation of the DMRG method (with a con- 
stant number of particles)^, and by functional renormal- 
ization group studies^. In these works it was argued 
that for the system sizes treated, the asymptotic regime 
in which the x^^ behavior is predicted is not yet reached. 
In Ref.[l3 it was shown that using the fRG method, which 
is argued to be as accurate as the DMRG method, even 
L of the order of 10^ is not sufficient to obtain the values 
of g of Eq. (O for I/t > 1. 

We now turn on the disorder by taking W ^ 0. In this 
case the results of ANj are averaged over 100 different 
realizations of disorder. In Fig. [51 the averaged particle 
density for = 0.1 is shown and compared to the W = 
case for various interaction strengths. As can be seen, for 
small values of the interaction the effect of disorder is very 
weak, while for large values of /, the FO decays faster in 
the presence of disorder. Zooming into these curves, it 
can be shown that the effect of disorder is to multiply the 
clean FO decay by an exponential factor e~^/^, where ^ 
is a characteristic decay length. 

For each strength of the interaction, one can rescale 
the disordered W curves, to the clean W — one 
by simply multiplying it by e^^^, using ^ as a fitting pa- 
rameter. As can be seen in Fig. [SI by using this rescaling 
method, the averaged disordered data collapses on the 
curves of the clean sample. 

The dependence of the decay length ^ on the interac- 
tion strength / is shown in the inset of Fig. [51 We shall 
now show that this quantity ^ is effectively the mobility 
localization length. 

The effect of disorder in the continuum limit can be 
divided to the forward and backward scattering terms. 
Whereas the backward scattering term is related to the 
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FIG. 5: The decay of FO for / = 0.0, 0.5, 1.0 and 1.5 (bottom 
to top) with L = 500 and eo = 10. The symbols are for 
W = 0.1 and the lines are for the clean {W = Q) case. The 
disorder effect becomes significant for large values of / where 
the localization length is small. The average was done over 
100 realizations. 

conductance and thus to the locahzation length of the 
electrons, forward scattering processes contribute only 
to the decay length of the FO, but not to localization. 
Thus, at first sight ^ does not necessarily correspond to 
the localization length. Nevertheless, in this case one 
can argue that the contribution of the forward scattering 
process to ^ is small and therefore ^ is a good measure 
of the localization length. 

Using standard bosonization technique it can be shown 
that the forward scattering processes result in the follow- 
ing term in the Hamiltonian: 

Hfs ^ - j dxi]{x) ^\7<P, (8) 

where ip is the TLL field which is related to the density 
operator by p{x) ^ —-^ '\7 'Pi^) ^-nd 'ri{x) is the q ^ com- 
ponent of the random potential. Since the TLL Hamil- 
tonian (u being the velocity) 

Htll^^ fdx[g{v^r + kvdn (9) 

depends on (f) only through {\/(t){x))'^ , it is easy to show 
that by a redefinition of the field = ~ f dyri{y) 
one can incorporate the Hfs term inside Htll and get 
a similar form of Hamiltonian. Therefore, the forward 
scattering term is not expected to change the physics of 
the system. 

Nevertheless, it was shown that this redefinition of the 
field has an effect on the correlation functions^. This 
results in a decay of the density-density correlation func- 
tion, which is, practically, the quantity we measure, and 
this decay is not related to the conductance. It is an ex- 
ponential decay of the form e~^/', where I = 2Sj(f)^' 
and Df is the forward scattering strength of the disorder 
(defined in the non-interacting case). 



For the decay described by the characteristic length 
one can find, using the Bethe Ansatz solution, the factor 
u/g for each value of /. It is easy to show that u/g, and 
thus I, are monotonically increasing functions of /, as 
opposed to the FO decay length (see Fig. [6]in the inset). 

Moreover, one can estimate I quantitatively for the sys- 
tem we deal with. The factor u/g found from Bethe 
Ansatz solution ranges from u/g = 2 for I — Q to u/g ^ 
4.5 for / = 1.5<. Denoting the amplitude of the disorder 
correlation function by D, i.e. {V{x)V[x')) = D5{x — x'), 
one finds that Df and Di, (the forward and backward 
scattering disorder strengths) are of the same order of 
magnitude as D. For non-interacting spinless electrons 
m a one dimensional lattice^^ l/Dt lOO/W^. Substi- 
tuting W = 0.1, one gets I of the order of 10^, which is 
much longer than the observed decay length. 

We thus conclude that the backward scattering pro- 
cesses are much more significant in the model treated, 
thus ^ is a very good approximation to the localiza- 
tion length, and its interaction dependence should be 
described by Eq. ([T]). 

Using the prediction of Eq. ([T]) with the value of the 
disorder we employ along this paper (order of 10~^), and 
recalling that without interactions ~ 100 /W^, the lo- 
calization length should range between ^(/ = 0) ^ 10'', 
which is much larger than the lattice sizes we considered, 
and thus almost doesn't infiuence the electron density, 
to (^(/ = 2) ^ 10^, in which the disorder effect should 
indeed be much more dominant, in agreement with the 
qualitative results presented in Fig. [5] 

The quantitative data (Fig. [S] (inset)) fits the theoret- 
ical predictions of Eq. ([1]) for not too weak interactions. 
For weak interactions (J ^ 0.5) no such fit was found, 
which is however expected, since for this regime the the- 
oretical localization length is much larger than the wire 
length. The fact that the best fit to Eq. ([T]) was for 
^0 ^ 7000 (not the expected Co ~ 10000) can be at- 
tributed to the same reason, as well as to the neglected 
forward scattering term which is stronger for weak /. 
We also note that the exact choice of the wire slices over 
which the fit is done, can change slightly the values of 
This, however, does not change the qualitative results, 
showing monotonic decrease of ^ as a function of the in- 
teraction strength. 

To summarize, the effect of disorder on the FO decay 
in the Anderson regime can be described by an extra 
exponential decay of the FO, which depends on the lo- 
calization length, of the form 

< AiV, Ai-iyr^ expi-j/ag)), (10) 

where the localization length ^(g) decreases monotoni- 
cally as the interactions increase. 

IV. DISCUSSION AND CONCLUSIONS 

In conclusion, we have shown that the FO envelope, 
in the TLL phase and the resulting disordered Ander- 
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FIG. 6: The rescaled decay of FO for the W — 0.1 curves over 
the clean curves with different interaction strengths. Inset: 
the locahzation length found by the best fit for each value of 
/ (symbols) . The line corresponds to the theory prediction of 
Eq. © with 5o = TOOO. 



a very clear FO, even for a specific sample. 



APPENDIX: FRIEDEL OSCILLATIONS IN THE 
ID TIGHT-BINDING MODEL 

In this appendix we calculate N{m), the density of 
electrons in site m, of a half filled one dimensional Tight- 
Binding lead, which is coupled to an impurity, in the 
asymptotic (to ^1) limit. The system is described by 
the Hamiltonian 



H = eocjco - V{cqCi + h.c.) 

L-l 



(A.l) 



N{m) can be calculated using the retarded Green func- 
tion of an electron in the m'th site G^{Lu;m,m)^, and 
the relation (for a half filled band) 



son insulator phase, is affected by both the interaction 
strength, and the disorder strength. Interactions actu- 
ally enhances the effect of FO since it drops with a weaker 
power law while disorder decreases the FO oscilla- 
tion since it adds an exponential factor to its decay form. 
We have shown that the length scale for this exponential 
decay is a good approximation to the mobility localiza- 
tion length, since it is weakly influenced by forward scat- 
tering processes for weak disorder. Thus we established a 
convenient way to evaluate the dependence of the local- 
ization length on disorder and interaction using only the 
ground state properties of the system. Qualitatively, the 
localization length as a function of interaction for a given 
weak disorder always decreases. As long as the localiza- 
tion length is not much longer than the wire's length the 
localization length behavior is quantitatively described 
by the renormalization group results^. Wc have also an- 
alytically described the dependence of the FO amplitude 
on the impurity strength. 

Finally we remark to the experimental relevance of this 
work. The theoretical treatment of disorder usually in- 
volves statistics over an ensemble of many samples which 
is usually hard to obtain experimentally. Furthermore, in 
the case we deal here, a measurement of FO on a disor- 
dered sample seems at first sight daunting. However, the 
simple method we suggest in order to deal with the dis- 
order, is in principle experimentally feasible, and solves 
these two difficulties. 

Once a technical method for measuring the electron 
density is established, it should be used twice for each 
sample, before and after the coupling of the wire to the 
dot. In principle, by using a gate it should be possible to 
eliminate the coupling between the dot and the wire. Our 
results, as can be seen in Fig. [51 which presents typical 
results for a particular realization, point out that the 
difference between these two measurements should show 



N{m) 



TT 



G^{uj] TO, rnjduj, 



(A.2) 



2t 



where we are possibly neglecting bound states with en- 
ergy lower than —2t, which give exponentially small con- 
tributions for large m. The Green function itself is de- 
termined by 

G^{uj;m,m) = G^{uj;m,m) + (A.3) 
G^{uj;m,l) ■ V ■ G(w;0,0) • V ■ G^(w;1,to). 

In this expression G^{Lo\m,l) is the bare (i.e., with- 
out dot) lead Green function, while G((jj;0,0) = {uo — 
eo — 2](w))~^ is the dot's Green function, where S(aj) = 
■!-(§- - is the self energy of the dot^^. The 

first term in the RHS of the equation simply gives the 
constant n — 1/2 occupation in the absence of the dot. 

Substituting the known wave functions and energies of 
the tight-binding Hamiltonian one finds 



1 



E 

/c>0 



cos(fca(TO — I)) 



(A.4) 

cos(fca(TO -I- /)) 



2t cos(fca) 



where k — J- fife, for integer n^. Transforming to integra- 
tion over unit circle in the complex plane leads to 



G«(c^;to,1) 



UJ 

2t 



^Wi-(-)^J 



(A.5) 



Combining Eqs. (|A.2|) . (|A.3|) and (|A.5|) . one can get. 



AA^(m) = iV(TO)-f/2 = 



(A.6) 



y2 nO 



LU - eo 



t \2t 
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and by substituting uj = —2tcos9, we find 



Y2 ,7^/2 

AN(m) = — / de- 



sin(6')e 



i2m0 



l-Tr/2 2t cos{9) + eo - 
One now defines z = e^*^ in order to get 



AN{m) 



z + eoz/t+l-V^/t'- 



(A.7) 



.,(A., 



where the integration is over the right half of the unit cir- 
cle, between the points ±1 on the imaginary axis (contour 
A in Fig.IZl). 



FIG. 7: The integration contours A and B-C-D which connect 
the points [0, —1] and [0, 1]. 



Next we deform our contour to the contour B-C-D in 
Fig. [7l In doing so we neglect the contribution of poles 
which may occur inside the closed line A-B-C-D. These 
represent states bound at the impurity, and as we have 
mentioned above, contribute exponentially small terms 
for large m. The integration in parts B and D is done 
by defining z = ±ix, respectively, x G [l,oo), while the 
contribution of the semicircle C vanishes as its radius 
goes to infinity. Therefore we get 



AN{m) = 
(-1)" 3 



(A.9) 



+ 1) 



TTf^^ ^' ^ ./i a;2+ieo.T/i-l + yVi2' 



-2m-l 



dx. 



For m ^ 1 the term x^2"i-i varies much faster than 
the other terms, and the rest of the integrand can be 
evaluated at x ~ 1 to give yrTp^p^j^- One thus gets the 
final form 



AN{m) 



nm eot' 



(A.IO) 
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